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Abstract 

We introduce a class of cellular automata associated with crystals 
of irreducible finite dimensional representations of quantum affine alge- 
bras UL(Q n ). They have solitons labeled by crystals of the smaller al- 
gebra t/q(0n-i). We prove stable propagation of one soliton for g n — 
42-1,4^ B£\CP, and D%l v For g n = C ( n 1} , we also prove that 
the scattering matrices of two solitons coincide with the combinatorial R 
matrices of t/^C^i ^-crystals. 



1 Introduction 



Cellular automata are the dynamical systems in which the dependent variables 
assigned to a space lattice take discrete values and evolve under a certain rule. 
They exhibit rich behavior, which have been widely investigated in physics, 
chemistry, biology and computer sciences Jwfl . When the space lattice is one 
dimensional, there are several examples known as the soliton cellular automata 
[|P|, [Pf], iPASj |PST| , |T|, 0. They possess analogous features to the solitons in 
integrable non-linear partial differential equations. For example, some patterns 
propagate with fixed velocity and they undergo collisions retaining their identity 
and only changing their phases. 

There is a notable progress recently in understanding the integrable structure 
in the soliton cellular automata. In the papers |TTMS| , [MSTTTj , [TNS| it was 
shown that a class of soliton cellular automata can be derived from the known 
soliton equations such as Lotka-Volterra and Toda equations through a limiting 
procedure called ultra-discretization. The method enables one to construct the 
explicit solutions and the conserved quantities of the former from that of the 
latter. A key in the ultra-discretization is the identities: (a, b <E R) 



lim elog(e< 
lim elog(e"< 



+ e«) = max(a, 6), 
b. 



x e< 



a 



'Institute of Physics, University of Tokyo, Komaba, Tokyo 153-8902, Japan 
^Department of Mathematics and Physics, National Defense Academy, Yokosuka 239-8686, 
Japan 



In a sense they change + into max and x into +. This is a transformation of 
the continuous operations into piecewise linear ones preserving the distributive 
law: 

(A + B) x C = (A x C) + (B x C) -> max(a, b) + c = max(a + c,b + c). 



The non-uniqueness of the distributive structure is noted by Schiitzenberger in 
combinatorics, where the procedure corresponding to the inverse of the ultra- 
discretization is called 'tropical variable change' [jKj. 

There is yet further intriguing aspect in the soliton cellular automata (called 
'box and ball systems') in JtJ, TNS|. There the scattering of two solitons is de- 
scribed by the rule which turns out to be identical with the U' q {An^) combinato- 
rial R matrix [NY| from the crystal base theory. The latter has an origin in the 
quantum afhne algebras at q — 0, where the representation theory is piecewise 
linear in a certain sense. 

Motivated by these observations we formulate in this paper and | HHIKTT ] a 
class of cellular automata directly in terms of crystals and link the subject to the 
1 + 1 dimension al qu antum integrable systems. The theory of crystals is invented 
by Kashiwara | Kas as a representation theory of the quantized Kac- Moody 
algebras at q — 0. It is a powerful tool that reduces many essential problems into 
combinatorial questions on the associated crystals. Irreducible decomposition of 
tensor product s and t he R obinson-Schensted-Knuth correspondence are typical 
such problems | DJM , |KN|, [n| . By connecting the classical and afhne crystals, it 

N2j 



also explains KMN1 



KMN2| the appearance of the afhne Lie algebra characters 



[DJKMO| in Baxter's corner transfer matrix method in solvable lattice models 

Here we shall introduce a cellular automaton associated with crystals of 
irreducible finite dimensional representations of quantum affine algebras U' q (Q n ). 
The basic idea is to regard the time evolution in the automaton as the action 
of a row-to-row transfer matrix of integrable U' q (Q n ) vertex models at q = 0. 
The essential point is to consider the tensor product of crystals not around the 
'anti- ferromagnetic vacuum' as in pCMNl] , |KMN2| ], but rather in the vicinity of 
the 'ferromagnetic vacuum'. 

Let B be a classical crystal of irreducible finite dimensional representations 
of the quantum affine algebra U' q (Q n ). It is a finite set having a weight decom- 
position and equipped with the maps e.;, fi : B — > B U {0} and <pi : B — > Z>o 
(i 6 {0,1,... , n}) satisfying certain axioms, (cf. Definition 2.1 in |KKM|.) For 
two crystals B and B' the tensor products B' ®B and B®B' are again crystals 
which are canonically isomorphic. The isomorphism B' <E> B — > B <S> B' is called 
the combinatorial R matrix Suppose that there are special elements denoted 

1 More precisely, the isom orphism combined with the data on the energy function is called 
the combinatorial R matrix [KMN1|. 
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by 1 G B and 6 B' with the properties 



(I) m ® 1 ~ 1 ® 

(II) For any u £ B' there exists fc € Z>o such that 

k 



u® 1 (g) • • ■ <g> 1 ~ b\ <g> ■ ■ ■ <8> 6 fc <g> mi,, (heB). 

We take dynamical variables of our automaton from the crystal i? and regard 
their array ... ,b t _ ll b t 0l b\ 1 . . . at time t as an element of the tensor product of 
crystals 

■■■ ®b t _ l ®b t ®b\® ® B ® B ® B ® ■■■ , 

where we assume the boundary condition &*■ = 1 € B for \ j\ 3> 1. See Section |2.2| 
for a precise treatment. The 'ferromagnetic' state V&* = 1 is understood as the 
vacuum of the automaton. The infinite tensor product with such a boundary 
condition does not admit a crystal structure. Nevertheless one can make sense of 
the construction below thanks to the properties (I) and (II) . The time evolution 
is induced by sending u§ from left to right via the repeated application of the 
combinatorial R matrix as 

B'®(---<8 B®B®B®- ■ ■ ) ~(- • -® B® B® B®---)®B' 
Ui,<8>(- ■ ■®b t _ 1 ®b t ®b\®- ••)-(•• •®fr*_ + i 1 ® & o +1 ® & i + V • ■ )®u^, 

which is well-defined as long as the above properties and the boundary condi- 
tions are fulfilled. In the language of the quantum inverse scattering method 



[3TF, KS|, this is the action of the q — row-to-row transfer matrix whose aux- 
iliary and quantum spaces are labeled by B' and ■■■®B®B<E>---, respectively. 
Note that the transfer matrix has effectively reduced to the (u^, W[,)-component 
of the monodromy matrix since its action is considered under the ferromag- 
netic boundary condition. The fundamental case g„ = will be studied in a 



more general setting in [HHIKTT]. In this paper we concentrate on the other 
non-exceptional series 

a - A {2) A (2) c« nW d (2) 



with the following choice of crystals: 



B = B 1 3l, B' = B^ 3 u^. 
Here B x is the crystal associated with the vector representation of the classical 

(2) (2) 

subalgebra of g n except for A\^ and D ^ Z^. Their cardinalities are %B — 2n, 2n+ 
l,2n + l,2n, In and In + 2, respectively. The element 1 e B\ is the highest 
weight one[j. To explain B^ and u^, recall the coherent family {Bi \ I 6 Z>i} of 

2 For (g n ,Bi) treated in this paper, a parallel construction seems possible also with the 
choice 1 = lowest weight element. 
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the perfect crystals obtained in [KKM|. It contains the B\ as its first member 
[]. The Bi with higher I corresponds to an Z-fold symmetric fusion of B\ . Then 
B^ in question is an infinite set corresponding to a certain I — > oo limit of Bi 
and is its highest weight element Q We shall call the resulting dynamical 
system U q (Q n ) automaton. They are essentially solvable trigonometric vertex 
models at q = in the vicinity of the ferromagnetic vacuum. A peculiarity here 
is the extreme anisotropy with respect to the relevant fusion degrees; B\ is the 
simplest one, while B' = B^ corresponds to an infinite fusion []. 

Once the automata are constructed the first question will be if they are 
solitonic. We prove a theorem that 

• the U q (g n ) automaton has the patterns labeled by the crystals {B{\ of the 
algebra U' q (Q n -i) that propagate stably with velocity I. 

Computer experiments indicate that they indeed behave like solitons. For in- 
stance, the initially separated patterns labeled by the L^(g„_i)-crystal elements 
b £ Bi and c € Bk (I > k) undergo a scattering into two patterns labeled again 
by some d S and b' e B{. Let S : Bi®B^ — > Bk®Bi be the two-body scatter- 
ing matrix of such collisions, namely, S(b®c) = d<g>b' . Let 1Z : Bi®Bk — > B^Bi 
denote the combinatorial R matrix of U q (g n -i). Then we prove 

• s = n 

for Q n — CrP and conjecture it for all the other §„. Similarly the scattering of 
multi-solitons labeled by B^, . . . , Bi N (?!>•••> In) is given by the isomor- 
phism B\ x % ■ ■ ■ <g) Bi N ~ Bi N Cg> • • • <g> Bi x experimentally. Thus the solitonic 
nature is guaranteed by the Yang-Baxter equation obeyed by S = 1Z. A precise 
formulation of these claims is done through an injection 



l^(j„_i)-crystal B t -» (l^(g„)-crystal B 1 



which will be described in Section ^. 

Admitting that they are soliton cellular automata, the second question is if 
there exist classical integrable equations governing them, possibly via the ultra- 

(2) 

discretization. Here we only confirm this for A 2 case by relating the associated 
automaton to the known example TS|. This observation is due to HI . 
The layout of the paper is as follows. In Section @ we first explain our 

("2") 

construction of the U q (Q n ) automata concretely along the Q n = A 2 example. 
It is valid for any U' q (Q n ) and any finite crystals having the properties (I) and 
(II). In Section |, we formulate the theorem and the conjecture for g n precisely. 
We sketch a proof of S = 1Z for case. In principle the idea used in the 
proof can also be used for the other g„. We will specify B^ as an infinite set 
with the actions e,, : B^ — > Bh U {0} but without the maps £j, ipi. In Section 



-,(D 



the family in [ KKM 



does not contain B 

and boo in [|KKM 



See [ HKKOT 



3 For C£ 

4 They are different from the limits B, 
5 To take B' = Bi with finite I is an interesting generalization. Sec [HHIKTT 

case. 



for A 



(i) 
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|] concluding remarks are given. Appendix [A] is devoted to an explanation of 
what is meant by 'B^ ® B% ~ B\ ® B^\ which shows up when the infinite set 
is substituted into the finite crystal B' . This is actually abuse of notation 
meaning an invertible map R' : B^ x B\ — > B\ x B^ between the sets. We state 
a conjecture on a stability of the combinatorial R matrix Bi ® B\ ~ B\ % Bi 
when I gets large, which ensures the well-definedness of the map R' . It assures 
that we may regard B^ as a finite crystal Bi with a sufficiently large / to define 
our automata. 



Our construction here and in | HHIKTT is a crystal interpretation of the 



L-operator approach | HIK | for a Q n = case. The U'(An) automaton in 

this sense coincides with the ones in JTS| [i], |TNS[ ] . As in the C$p case in this 
paper, the properties stated in the above can actually be proved by means of 
the crystal theory. The detail will appear elsewhere along with the results on a 



more general choice of the crystals B and B' | HHIKTT | 



Acknowledgements The authors thank K. Hikami, R. Inoue, T. Nakashima, 
M. Okado, T. Takebe, T. Tokihiro, Z. Tsuboi and Y. Yamada for valuable 
discussions. 



(2) 

2 A 2 example 



(2) 

Let us explain our automata concretely along the case g n = A 2 . This simple 
example is helpful to gain the idea for the general g n case treated in the next 
section. 

(2) 

As a peculiarity in the ran k I situation, the U' q (A 2 ) automaton turns out 
to be an 'even time sector' of iTSf . 



2.1 A^ ] crystals 

For I E Z>i set 

Bi = {{x,y) | x,y E Z> , x + y < I}. (2.1) 
The action of the Kashiwara operators ej, fi(i = 0, 1) : Bi — ► Bi U {0} are given 

by 

eo{x,y) = < e 1 {x,y) = (x + l,y-l), 

[{x,y + lj nx<y, 

it \ )(x + l,y) Hx>y, - 

fo{x,y) = <, 1S .. fi{x,y) = {x-l,y+l). 

{(x,y-l) iix<y, 

In the above, the right hand sides are to be understood as if they are not in 
B x . Setting £,(6) = uiax k {^b ^ | k > 0} and <pi{b) = uia^ k {ffb ^ | k > 0}, 
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one has 



£oO) = I - x - y + 2(x - y)+, e 1 (b) = y, 
<Po(b) = 1 - x - V + 2(2/ - x )+, ( Pi(b)=x, 

for b — (x,y) S Here the symbol (•)+ stands for 

(x) + = max(0, x). 

(2) J L 

These results are obtained by extrapolating the result KKM | to n = 1. 
For ! = 1 we use a simpler notation as 



B = B U 1 = (1,0), 2 = (0,0), 3 = (0,1). 



(2.2) 



Given t wo c rystals B and one can form another crystal (tensor product) 
B ® B' \ Kas | . The crystal Bi ® B\ is connected and so is B\ ® B\. Calculating 
the map Bi® B\ — > By® Bi commuting with e< and one has 

Proposition 2.1. The combinatorial R matrix Bi ® B ~ B ® Bi is given by 



(x,y) ® 1 



(x,y)®2 



(x,y) ® 3 



1® (2,0) if {x,y) = (1,0), 

3®(x + l,y-l) ifx + y = l,y>l, 

2®(x + l,y) ifx + y = l-l, 

1 (8 (a; + 1, j/ + 1) otherwise, 

18(2-1,0) i/(x,y) = (Z,0), 
3®(x,y-l) ifx + y = l,y>l, 
otherwise, 

if (x,y) = (0,1), 
if (x,y) = (0,1-1), 
if(x,y) = (1,0), 
if2<x<l, y = 0, 
ifx = 0, 0<y<l- 
otherwise. 



2® (x,y) 

'3® (0,1) 
2® (0,0 
1®(0,1) 
1® (a: -2,y) 
1 ® (x, y + 2) 
3 ® (x — 1, y - 



1 



In Section 2.2 we shall use formal / — > oo limits of and the combinatorial 
i? matrix Bi ® B ~ B ® Bi. In the present case the prescription is to simply 
shift the coordinate (x,y) to (x — l,y) and to consider 

= {(z, y) | x e Z< , y e Z> , x + y < 0}, 

without specifying a crystal structure. The map R' ; B^ ® B ~ B ® in the 
sense of Appendix |A] is deduced from Proposition 2J by concentrating on those 



G 



(x,y) in the vicinity of (0,0). Thus it reads 

•l®(0,0) if {x,y) = (0,0), 

{rin , 3^ + 1,,-!) ifx + y = 0y>l, 

2 <g) (a; + 1, y) \ix + y = —l, 

1 <8> (a; + 1, y + 1) otherwise, 

'l®(-l,0) if (x,y) = (0,0), 
(.<•.//) • <( :5 i,r. // 1 i ifx + y = 0,y>l, (2.4) 

2 ® (a;, y) otherwise, 

MgjJ 10 !*-^ , * V = > (2.5) 
13(8) (a; — l,y — 1) otherwise. 

To depict this in a figure we put fc b , beside the arrow (x,y) — > (x',y') to signify 
the relation 

i?' : (x,y) (8) 6 6' (g> (a/,y')' 



We call 6 and 6' the upper index and the lower index, respectively. Now (2.3)- 
( [2.5] ) are summarized in the semi-infinite triangle in Figure [l| 

2.2 Cellular automaton 

By applying B^® B ~ B ® successively one has 



B^B <g>---®.B~B(gi---<8 fl(8»Bi, 
u<8>&i <g> • • • <g> <8> • • • ® &^<8 u' (2.6) 

for any L G Z>j , where i = — [ 4], j = i + L—l f\ Denote the element (0, 0) G B^ 
by ut,. The map R' : B^ ® B ~ B ® B^ (^3|)-(p^) has the properties (I) and 
(II) in Section 0. Set 

7? = {fe : Z -> B | bj = 1 G B for \j\ > 1}. 

We shall regard V as a subset of the tensor product which is formally infinite 
in both directions, i.e., 

V = {■ ■ ■ <g> 6_i <g> b <g> 6i <g> • • • e • • • <g) B <g> B ig) B ® • • • | b 6 = 1 for \ j\ » 1}. 

(2.7) 

In the latter picture one should distinguish the elements even though they are 
the same under translations. For example, 

•••l(8l(g>l<g>2(83<g>l<g)l<g>l<g>--' , 
•••l(8l<g>l<g>l<8)2<g)3<g)l<g>l<g>--' 



6 The symbol [x] denotes the largest integer not exceeding x. 
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Figure 1: 

The Semi-infinite triangle representing R' . There are 6 different patterns de- 
picted by circles, squares and diagonal squares which are filled or empty. 



are distinct elements in V. The set V ( |2.7| ) is not equipped with a crystal struc- 
ture. Nevertheless the properties (I) and (II) enable us to define an invertible 



map T : V — > V that formally corresponds to an L — > oo limit of (2.6). To 
describe it precisely, note that any element in V has the form 

p = ---®l®l<g>6i<8>---i8)&j-<8)l<8>l<8>--- (i<j), (2.8) 

where 6^,6^+1,... ,6j € B. Owing to the properties (I) and (II) there exists 
fco G Z>o such that 



U[, ® 6i <g> • • • <g> 6j- <g> 1 ® • • • <g> 1 ~ b'f ® • • • ® bj <8> 6j- +1 <g> • • • (8 6^ +fc <8> 

for all fc > fc - Then T(p) is defined by 

T(p) = • • • <g> 1 <g> 1 <8> b'i <g> • • • <g> ^ (g) ® • • • <g> ^ +fe (g) 1 <8> 1 ® • • • , 

which is fc- independent as long as k > fco- The inverse T _1 can be described 
similarly. 

The map T plays the role of the 'time evolution' operator. It is a q = 
analogue of the row-to-row transfer matrix of a solvable lattice model in the 
vicinity of the ferromagnetic vacuum. 

Given p € V in fl2.8| ) define u m G for all m G Z by the recursion relation 
and the boundary condition 

Ufn-l <8> &m — Kn ® w m for all 771 G Z, 
7i m = lt|} for 777 < 7 — 1. 

Plainly, «(,(£>(•• •<8& m -i<8>6 m (g>& m+ i<8>- ••) — •" ' ® ® Um ® kro+i ® ' ' ' • 

Due to the properties (I) and (II) the sequence u m , u m +i, - - . tends to = 
(0,0) G -Btj. In this way any element p G P specifies a trajectory {um}m=-oo m 
the semi- infinite triangle (Figure [IJ) that starts at the origin (0, 0) and returns to 
it finally. This picture is useful in calculating T(p). Namely, the trajectory is de- 
termined by following the arrows with the upper indices . . . , 6 m -i, b m , 6 m +i, . . . 
appearing in p = ■ ■ ■ <g> 6 m _i (g> 6 m ® 6 m +i ® ■ • • ■ Then T(p) is constructed by 
tracing their lower indices as T{p) = ■ ■ ■ <g) ® 6^ <8> ® • • • • 

For p = • • • <g> bj <8 6j+i and i G Z define 6* G B by T*(p) = 

• • • <8> bj <8> b* +1 ® • • • . Then the time evolution of the cellular automaton is 
displayed with the arrays 



6° 


to 


b° 


1 


b° 


b\ 


bl 


b 1 


2 


b 1 


1 


% 


b\ 


b\ 


b 2 


2 


b 2 


1 


bl 


b\ 


bl 



Let us present a few examples. 
Example 2.2. 
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t=0 

t=l 

t=2 
t=3 



• 11233311111111111111111111111- • • 

• 11111111123331111111111111111- • • 

• 11111111111111112333111111111- •• 

• 11111111111111111111111233311- • • 



Example 2.3. 



t=0 
t=l 
t=2 
t=3 



11333311111111111111111111111111- • • 
11111111113333111111111111111111- - ■ 
11111111111111111133331111111111- • • 
11111111111111111111111111333311- • • 



Example 2.4. 



t=0 : •••1133311111123111111111111111111111 
t=l : •••1111111133311123111111111111111111 
t=2 : •••1111111111111133123311111111111111 
t=3 : •••1111111111111111123111133311111111 
t=4 : •••1111111111111111111123111111133311 



Example 2.5. 

t=0 : - - - 11233111331111111111111211111111111111111111111- ■ ■ 
t=l : - - - 11111112331133111111111121111111111111111111111- ■ ■ 
t=2 : ■ - - 11111111111133112331111112111111111111111111111- ■ ■ 
t=3 : ■ - - 11111111111111113311123311211111111111111111111- ■ ■ 
t=4 : ■ - - 11111111111111111111331111223311111111111111111- ■ ■ 
t=5 : ■ - - 11111111111111111111111133121111233111111111111- ■ ■ 
t=6 : ■ - - 11111111111111111111111111121331111112331111111- ■ ■ 
t=7 : ■ - - 11111111111111111111111111112111133111111123311- ■ ■ 

Example 2.6. 

t=0 : ■ - - 11233111111133111121111111111111111111111111111- ■ ■ 
t=l : ■ - - 11111112331111113312111111111111111111111111111- ■ ■ 
t=2 : ■ - - 11111111111123311112133111111111111111111111111- ■ ■ 
t=3 : ■ - - 11111111111111111233211113311111111111111111111- ■ ■ 
t=4 : - - • 11111111111111111111211233111331111111111111111- - • 
t=5 : ■ - - 11111111111111111111121111112331133111111111111- ■ ■ 
t=6 : ■ - - 11111111111111111111112111111111133112331111111- ■ ■ 
t=7 : ■ - - 11111111111111111111111211111111111113311123311- ■ ■ 

The last two show the independence of the order of collisions. These 
pies suggest that the following patterns are stable (Q £ Z>i, R = 0, 1): 



Q 



■•■®2<8>3<g>---<g>3<8> 



R=l, 



Q 



•••<g>3<g>3<g)---<g)3<S) 



R = 0. 
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The both patterns should not be followed by 3. The former pattern can be 
preceded by any element in B while the latter should only be preceded by 1. Q 
is the size of the soliton and R is the number of occurrences of 2 in its front. 
They move to the right with the velocity 2Q — R when separated sufficiently. 

(2) n — I 

These features are consistent with g n = A 2 case of Theorem 3.1, See also 



Section 3.2. 

In fact the U q (A\ ) automaton described above can be interpreted [HI] as 
an 'even time sector' of the automaton in [TS]. Replace the array of {1, 2, 3} by 
that of {1, 2} with double length via the rule 1 — > 11, 2 — > 12 and 3 — > 22. In the 



resulting array, play the 'box and ball game' as in [TS|. Namely, we regard the 
array as a sequence of cells which contains a ball or not according to the array 
variable is 2 or 1, respectively. In each time step, we move each ball once to the 
nearest right empty box starting from the leftmost ball. Then the 2 time steps 

( 2} 

in the box and ball system yield the 1 time step in our U q (A 2 ) automaton. 

In terms of crystals, this can be explained as follows. First, the box and ball 
game in [TS] is known [ HHIKTT | to be equivalent to the U' q (A^) automaton. 
Let 



Bk = {mi . . . m k | mi £ {1, 2}, mi < ■ • • < m k } 

denote the U q (A^)-cijstal corresponding to the fc-fold symmetric tensor rep- 
resentation. (We have omitted the frame of the usual semistandard tableaux.) 
Consider the maps and h\ defined by 



: 

(x,y) 

h x : B 1 
1 
2 
3 



k-2y 2y 



B k (fc>l) 

k-\-x-\-y —x — y 



-> 1...12. 

B l ®B l 
1(8)1 
I® 2 
2® 2. 



.201. ..12. 



Then for k large enough, we have the commutative diagram: 

h^hi 



B^ (giBi 

R' 



B k ® B k ® B x ® B x 



B, 



h 1 m^ B l ®B l ®B k ® B k , 



where the down arrow in the right column is the crystal isomorphism. This 
asserts that the square of the T in the U'^A^p) automaton coincides with the 
T of the U' q (A { 2 2) ) automaton. 
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3 U'(g n ) automaton 



3.1 Theorem and conjecture 

Let us proceed to the U'(g n ) automata for g„ = A, 



(2) 

2n-l 



(n > 3),4 2 n } (n 



*2n V"- > 

(n > S),^ 1 ' (n > 2),£)^ 1) (n > 4) and D^ 2 ^ (n > 2). Our aim here 
is to formulate the theorem and the conjecture stated in Section FA precisely. In 

(2) 

principle construction of the automata is the same as the A 2 case explained in 
Section 2.2. The time evolution operator T is constructed from the invertible 
map R' : g B\ — > B\ % B^ in the sense of Appendix [A]. However its analytic 
form like (2.3)— (2.5 ) is yet unknown for general g„ []. Thus we have generated 
the combinatorial R matrix R : B\ ® B\ ~ Bi eg) B\ directly by computer, and 
investigated the automata associated with R instead of R! for several large I. 



Consistently with Conjecture A.l, their behaviour becomes stable when I gets 



large, yielding our automata associated with R' . 



What we present in Section B.2 is the list of the data Bi,Bh,u^ and ii for 
each Q n . We change the notation slightly from Section |l| and representing 
elements in B — B\ with symbols inside a box or (f>. For example the special 
(highest weight) element 1 € B\ in the properties (I)— (II) is denoted by [T]. 
Consequently a state of the automata at each time is represented by an element 



V = {b : Z -> B | bj = |T| S B for \j\ > 1}, 
which is also interpreted as 

V = {■ ■ ■ (8) 6_i ® b <8> bx ® • • • G • • • ® B ® B ® B ® • • • | bj = |T| for \j\ > 1}. 



We let T : "P — > denote the time evolution operator as in Section 2.2. The 
essential data is the injection 

H ■ ^(§„_i)-crystal B l -» ([/^( fln )-crystal 

which will be utilized to label solitons in the U' q (Q n ) automata in terms of 
elements of the L^(g ra _i)-crystal Bi. 

First we claim stable propagation of the 1-soliton as 

Theorem 3.1. For any I 6 Z>i and b £ U' q (Q n -i)-crystal Bi, the time evolu- 
tion map T acts on the state 

- * * <8> [Tf <s» [T] <g» «/(&) <s> [T] <s> [T[ <s» • * - 

as the overall translation to the right by I slots. 

Thus I is the velocity of the soliton. Even when n is the minimal possible 
value for g n , Theorem 3.1 makes sense if i\ and B\ for 'J7'(0„_i)' are interpreted 



7 For C„ we have a concrete descript ion of t he combinatorial R matrix B; ~ ®Bi 

in terms of an insertion algorithm. See [ HKOT ] . 
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appropriately by an extrapolation of the data given in Section 3.2. For example, 

(2) I ■ ' ' 

the solitons for g n = A y 2 case argued in Section 2.2 agree with those coming 



12) 

from iT I'(A Q )-crystal B{ in such a sense. The proof is given in Section 3.3 



9 

For scattering of two solitons we have 

Theorem 3.2. Let g n = Cn (n > 3). For fixed positive integers I > k, let 
Bi®Bk 9 bi®02 — C2®c\ G Bk®Bi be an isomorphism under the combinatorial 
R matrix of the U' q (Q n -i)-crystals . Then for m 3> I, there exists to > such 
that for any t > to and some m' 



T* : ■• ■ <E>[j]<g)[7]<EUi(&i) ®[T] ® «fc(Z>2) <8> |T| <g> |T| O |T| 

^ ® |T| ® |T| ® |T| ® « fe (c 2 ) ®|T|® m ® ?z(ci) ®|T|®|T|® ••• . 

In other words, we have the combinatorial R matrix as the scattering matrix 
of the ultra-discrete solitons. Compared with «z(&i), 1^(02) is shifted to the right, 
but we do not concern the precise distance. A sketch of a proof of Theorem [j.2| 



will be given in Section 3.4 



In fact we have a conjecture on iV-soliton case for general g„. 

Conjecture 3.3. Let N 6 Z>2- Fix positive integers k\, . . . , k^ (k\ > ■ ■ ■ > 

kfj) and the elements bi € B^ (i = 1,... , N) of the U' q (Q n ~i)-crystals. Define 
Cj £ Bfci by b\ <g> • • • ® bjsr ~ cpj ® ■ ■ ■ ® c\ under the isomorphism B^ ® B}~ 2 ® 
■ ■ -<S>Bk N — Bk N <8> ■ • ■<S)Bk 2 ®Bk 1 ■ For mi, . . . , mjv-i 3> k%, there exists t > 
such that for any t > to, 



0z, , I \<&>rn 1 I iQ^mTv-l /T * I 1 



' ^ ■ • • ® |T|® tk N (c N ) ®[T]® mi ® • • • ®[T]® mjv 1 ®« fcl (ci)®|T|®--- , 

ZioZcfc for some m^, . . . , m' N _ 1 . 

In particular, ci, . . . , Cat do not depend on to^'s (i.e., the order of collisions) 
as long as mi 3> k\. Again we do not concern the precise distance between 
«fci(&i) and %k N (cN) in the abo ve. 

The rank n in Conjecture |3.3| should be taken greater than the minimal 
possible values, i.e., A^-i (n > 4),A 2 2 „ ) (n > 2) 1 B ( n 1) (n > 4),d 1) (n > 
3),£>1 1} (n > 5) and D^ ] +l (n > 3). We have checked the TV = 2,3 cases by 
computer with several examples for Q n = , , ylg 2 ^ , , , and 
the TV = 3 case for cf 5 . 

Let us present a few examples from the U'A C^) automaton. We use the 
notation that will be introduced in Section p^. T he dynamical variables are 
taken from the crystal B\ = {|~T|, [2], [~3~|, ["§], | 2 1, | 1 1}. In the following examples 
we drop the boxes for simplicity. 

Example 3.4. 
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t=0 : • • • 1122321111133111111111111111111111 • • • 

t=l : ■ ■ ■ 1111112232111331111111111111111111 ■ ■ ■ 

t=2 : • • • 1111111111223213311111111111111111 • • • 

t=3 : • • • 1111111111111123333311111111111111 • • • 

t=4 : • • • 1111111111111111133123331111111111 • • • 

t=5 : • • • 1111111111111111111331112333111111 • • • 

t=6 : • • • 1111111111111111111113311111233311 • • • 

Compare this with 

(1, 1, 0, 2) ® (0, 1, 1, 0) ~ (0, 2, 0, 0) ® (0, 1, 2, 1) 

under the isomorphism B±® B 2 — B 2 ® B4 of U' q {C 2 )-crystals. 

Example 3.5. 

t=0 : • • • 11122311111232111111111111111111111111 • • • 

t=l : • •■11111111223111232111111111111111111111 •■■ 

t=2 : • •■11111111111112311223211111111111111111 •■■ 

t=3 : • ••11111111111111111311122232111111111111 ••■ 

t=4 : • • • 11111111111111111111311111222321111111 • • • 

t=5 : ■ • • 11111111111111111111111311111112223211 • • • 

Compare this with 

(0, 1, 0, 2) ® (1, 0, 1, 1) ~ (0, 1, 0, 0) ® (1, 0, 1, 3) 

under the isomorphism B 5 ® B 3 ~ B 3 ® B 5 of {/^(C^^-crystals. 
Example 3.6. 

t=0 : • • • 1113333111111231111123111111111111111111111111111111 • • • 
t=l : • • • 1111111113333111123111231111111111111111111111111111 • • • 
t=2 : • • • 1111111111111113331112332311111111111111111111111111 • • • 
t=3 : • • • 1111111111111111111133113311222211111111111111111111 • • • 
t=4 : • • • 1111111111111111111111111133331111222211111111111111 • • • 
t=5 : • • • 1111111111111111111111111111113333111111222211111111 • • • 
t=6 : • • • 1111111111111111111111111111111111333311111111222211 • • • 

Compare this with 
(0, 2, 2, 0) ® (0, 1, 0, 1) ® (0, 0, 1, 1) ~ (0, 0, 0, 0) ® (0, 2, 2, 0) ® (1, 0, 0, 3) 

under the isomorphism B & ® B4 ® B 2 — B 2 ® B4 ® B 6 of {/^(C^^-crystals. 
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3.2 Data on crystals 

Let us present the data B t , B^u^ and %i for each a„. Bi and Bu will b e specified 
only as sets. The crystal structure of Bi is available in [ KKM |. (See [ HKKOT | 
for Cn^ case.) For all the g„, £ B^ is given by 

= (Vx;, x; = 0) 

in the notation employed below. As a set Bi can be embedded into B^ by 

9i ■ Bi — > B^ 

(xi,Xi)^ (x, - lSu,Xi). 

Obviously any u G B^ has the inverse image in each B[ if I is large enough. It 

is easy to see that the composition q _£} y ^ fiij^ ^ y |g| 9^ & y |q| 

is independent of I when Z is large enough. Here we understand that gi(0) = 0. 
Thus one can endow B^ with the actions of ii , fi defined by this with I sufficiently 
large. (However it will not be used in this paper.) 

n 

Bi = {(xi, . ..,x„,x n , ...,xi) £ Z 2n \xi,Xi > o, y^X x i +*») = 0: 

1=1 

n 

B^ = {(xi, . . . ,x n ,x n , ...,xi)e Z 2n \xi < 0, Xi > (i ^ 1), Xi > 0, ^{Xj + Xi) = 0}. 



For B\ we use a simpler notation 



\X\ , . . . , X n , X n , . . . , X\ J 



[7] if Xj = 1, others = 0, 
[T] if x* = 1, others = 0. 



For b = (xi, . . . ,£„_!, x„_i, . . . ,xi) in l^(A2„_ 3 )-crystal 



«i( fo )=|_2j <8> • • • <S> [fi| ®[nj ®---<g>|_2j 



a -4 (2) 



Bi = {(xi,. . . ,x„,x n , . . . , xi) e Z 2n | Xj,Xj > 0, ^(x^ + X*) < /}, 

n 

i?t, = {(xi, . . . , x n , x n , . . . , xi) E Z 2 ™ xi < 0, Xi > (i 7^ 1), Xi > 0, 2j(x* + Xj) < 0}. 
For B\ we use a simpler notation 

{ [7] if Xi = 1, others = 0, 
[T] if Xi = 1, others = 0, 
4> if Xi = 0, xi — for all i. 
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For b = (xi, . . . ,x n -i,x n -i, . . . ,xi) in l^(A2„_ 2 )-crystal Bi, we define s(b) = 
T:=i 1 (xi+*i),s'(b) = [(l-s(b))/2}. 

If I — s(b) is odd, 

= 0(g)jT| <8>\2j ® •••«>[§ ®---®[2j <S> [T] 

otherwise 

. . r^|(g>s'(f>) r^ilglXl I |®X„_i I I ®X„_1 I 10X1 I |®s'(i>) 

*i( fo )=|T| ®[2] ®0 (81 [2] o|T| . 

When n = 1, the above notation for £?i and ( |2.2[ ) are related by 1 = IT], 2 = <f> 
and 3 = \ The solitons and their velocity mentioned in the beginning of 
Section |2.2| agree with the n = 1 case here. One interprets s(b) = and 
s'(b) = [1/2] for b from '[/^(Ap 2 -*) '-crystal £>;. Then, under the identification 
R = (1 - (-l)')/2 and Q = i? + s'(6), the velocity is indeed I = 2Q - R. 
Qn = : 
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Bi = \ (x 1: . . . ,X„, SCO, 55n, ■ ■ ■ ,Xi) G Z 2 " X {0, 1} 
II, ■ ■ • ,X n ,X ,X n , . . . ,X X ) G Z 2 " X {0,1} 

For B\ we use a simpler notation 

ix \ , . . . , X n , Xq 7 X n , . . . , X\j 



xq — or 1, 2;^, a;i > 0, 1 

2^0 +J2i=l( x i = 1 J' 

xq = or 1, xi < 0, 
Xi>0(i^ 1), Si > 0, 

5=0 + J2i=l( X i +Xi) =0 



if x, = 1, others = 0, 
[T] if Xi = 1, others = 0. 



For b = (xi, . . . ,x„_i,x ,x„_i, . . . ,x x ) in [^(B^J-crystal S/, 



|4f2'i 



I |®x„_i I |Q9X I |Q9x„_i | 1 

®[0] ®0 ®---®[2] 



®xo I — i(8)x„_i 



a -r (1) 



= |(xi, . . . ,x„,x„, . . . ,xi) G Z 2n 

B% = |(xi,...,x„,x„,...,xi) G Z 2 " 
For B\ we use a simpler notation 

(xi , . . . , X n , X n , . . . , Xi ) 



Xi , Xi ^ 0, I 

«>E l "=i( a; i + ^) e/-2Z J' 

xi <0, Xi > 0(t^ 1), x 4 > 0, ] 

o>££=i(&<+5i)e2Z I 



pT| if Xi = 1, others = 0, 
[T] if Xi = 1, others = 0. 
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For b = (xi, . . . ,x n -i,x n -i, . . . ,xi) in [/^(C^-J-crystal Bi, we define s(b) = 
E^(xi+Xi),s'(b) = (l-s(b))/2. 

-\®s'(b) r-=-i8>xi |— rn®2„-i I — i<8>xn-i r7T\® X1 rrn^s'C)) 



|— n 59Xn-l ^ I 1 C 



Qn = : 



= j(xi, . . . ,x n ,x n , . . . ,xi) G Z 2 ™ 

= . . . ,x n ,x n , . . . ,xi) G Z 2 ™ 

For i?i we use a simpler notation 

(Xl , . . . , X n , X n , . . . , X\) 



x — or Xfi — 0, x% - x 2 0, 

x n = or x n — 0, Xi < 0, 
x 4 > 0(i # 1), > 0, 
E"=i(aJ» +^») = 



[T] if Xj = 1, others = 0, 
[T] if Xj = 1, others = 0. 



For b = (xi, . . .,x n -i,x n -i, . . . ,xi) in U'JD^l ^-crystal B u 



®x„_i I i®x„_i 

n 



)0 



(8>xi 



0n = D% : 



5; = |(xi, . . .,x n ,x ,x„, . . . ,Xi) e Z 2n X {0, 1} 

B\f = |(xi, . . . ,x„,x ,x„, . . . ,xi) G Z 2 ™ x {0, 1} 
For £?i we use a simpler notation 

Si 3 (xi, . . . ,x„,x ,x„, . . . ,xi) = < 



x = or 1, Xi, Xi > 0, I 
+ E"=i( :z: » +**) < 1 J ' 

xo = or 1, xi < 0, 
x, > 0(i ^ 1), Xi > 0, 
+ YJi=li x i + < 



[T] if Xi = 1, others = 0, 
pT| if Xi = 1, others = 0, 
</> if Xi = 0, Xi = for all i. 



For 6 = (xi, . . . , x„_i,Xo,x„_i, . . . ,xi) in J7g(Z)^)-crystal Bi, we define s(b) = 
YZ=i\xi+Xi)>s'(b) = [(l-s(b))/2}. 



If I — s(b) is odd, 
otherwise 

*i(6) = [T| <»■■ 



igii I 1 (8s' (6) 



9x„_i 1 1(81x0 1 i®a:n-i 1 iQ9Xi 1 iq 

®[o] ®[n] <8> ■ ■ ■ <8>[2] o|T| 

l®x 1 i(8x„_i I 1(8x1 I 1 81s' (6) 

1 1 ® •■•(81 [2] ®|T| 



I |8>x„_i 1—181x0 I iC 

)|_nj <8>|_0j ®[n] 
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3.3 Proof of Theorem O 



Our proof below uses the crystal structure of B\ given in jfiKKOTj for g n = C n 1} 
and in [KKM] for the other types. 



Lemma 3.7. For i G {1, 2, . . . , n — 1} ; the following diagram is commutative: 



U'Mn-i) -crystal B t 



U'(g n -x)- crystal Bi U {0} 



U q (g n )-crystal Bf 



C/'(fl„)-cn/staZ Bf U {0} 



XTie same relation holds between fi and fi+\. 

The Kashiwara operators of U q (g n -i) and l^(g„)-crystals should not be con- 
fused although we u se the same nota tion. The proof of the lemma is due to the 
explicit rules for e-i [KKM, HKKOTf and the embedding of [/^(§ n _ ^-crystal B\ 



into t/g(§„_i)-crystal Bf as J7q(g n _i)-crystals KN |. Here g„_i is the classical 

subalgebra of g„_i. According to g„ = A^, b£\ 

it is given by g„ = C n , C n , B n , C n , D n and B n respectively. 



and D { ? +1 



Proof of Theorem 3.1. 

a -a {2) ■ 

(2) 

Since B\ of [/^(A^^-crystal is isomorphic to B(lAi) as C/ 9 (C n _i)-crystals, 
for any b G Bi there exists a sequence i\,...,i p G {1, . . . , n — 1} such that 
bo = &i ■■■ei 1 b with 6 = (Z,0, ...,Q) G B\. From Lemma B/7, i;(&o) — 
e ip+ i ■ ■ •e il+ l^^(6) is valid. 
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Taking sufficiently large L and f = /| • • • fn—xfnfn—i '"fit one nas 



u L ®n(b) ®|T| 

e< p +i— e il+ i | 



[if ' ® u L 

p+i — e 4l+ i I 
[if «Ui(&o) ® 



(L,o,...,o)®^®0 0i 0Z ®^®(i,o,...,o) 



(L,0,...,0)®[|p'®[lf 



-L JX+2I 
e 11 



(0,...,0,£)®[lf 

'1 

(o > ...,o J L)®[T] 8 '®[If 

(L,O,...,O)®0 02i 



0®'®^®(i,O,...,O) 

[Tf% (o,i,o,...,o,i-o 

<{ 

\Tf 2l <E> (0,...,0,l,L-l) 

ei + l fo L+2 '{ 

g® 2i ®(L,o,...,o). 



Thus the isomorphism B L ® Bf 21 ^ Bf 21 ® B L sends u L <g> «/(&) <8> [1]®' to 
|Tf n(b) <g> ul. This proves Theorem 3J. 



Let 6j, mj =(m,0,...,0)eB| of ^(A^_ 2 ), where < m < I. First we consider 
b = b^ case. Taking sufficiently large L and f = • • • fn-if„fn-i ' ' ' /ai one 
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has 



«x <8> **(6o°) <8>IXf 



(L,o,---,o)®[2f l ®[if l 

4 

(£,0,---,0)®[lf'®[lf 
(L,0,...,0)®^®^ 

'1 

(L,o,...,o)®[lf' ®[Ip 

e l Jo n 

(i,o,...,o)®0 02/ 

Thus the isomorphism S L ® J3f 2Z 



[T] s, 0[l] s, »(i,o o) 

'I 

g 0Z ®^®(L,o,...,o) 

ef/ 2, /T+ 2! | 

[if 21 ® [L- 1,1,0,..., 0) 
f 



[I] ® (L -1,0,... ,0,1,0) 

et +! /o 2, /T +2 '| 

^ [lf%(L,0,...,0). 

Sf 2 ' ® B L sends u L ® t,(6^) ® [7]®' to 
[Tj 8 ® «i(&o^) ® ml, verifying Theorem |3.l| for b = b$\ Explicitly it reads 

(L, 0, . . . , 0) ® [if ® [if * [if ® Qif ® (L, 0, . . . , 0). (3.1) 

Let us proceed to arbitrary t e 5| case. We will attribute the proof to 
the b = fog™ 1 ' case for some < m < I and eventually to the bfp case shown 
above. Set k = [l=p], k' = [i=^±i] so that n(b ( ™ ] ) = ^ k '~ k ®[lf ®[lf m ® 

|Tf according to Section |3~2] . Since i?; of U' q (A^_ 2 )-crystal is isomorphic to 
B(/Ai) e B((l - l)Ai) © • • • © 5(0) as f7 g (C„_i)-crystals, for any b e Bi there 
exists a sequence ii, . . . ,% v G {1, . . . ,n — 1} such that 6q ' 
< m < I. From Lemma 3.7 one has 



ejj 6 for some 



ML ® ll(b) ®|jj 
e ip+ i— e il+ i I 

Ui ®^(6( m) )®0 0i+fc ' 

rk+k' rL + k+2k' 

e l Jo Ji I 

(Z,0,...,0)®[]L]®%[^®[Tf 



ML 



[if ®»i(&)®[ir 

[ip , ®» I ($ B) )®[ir*'®«£ 



S l Jo Ji I 



®[[f ®Q[f ®OM,...,0), 

where in the bottom we have used (3.1) and ui ® [T] — [T] ® ml. Thus the 
isomorphism B L ® Bf 21 ^Bf 21 ® B L sends u L ®«z(&)®[lf to [T] <S>Z ® n(b) ® u L , 



proving Theorem 3.1 
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n - R«- 
Q n — Dn . 

The proof is similar to the case g, 

a -r (1) - 

The proof is similar to the case g, 

fln =D [ n ] : 

The proof is similar to the case g, 
The proof is similar to the case g, 



41 x with f = /! • • ■ fUfn'fl-l ■■■II 

A^_ x withf = /l---/L2/^-i---4 

4 2 J withf = ^-..^_ 1 ^_ 1 .-./i 



3.4 Proof of Theorem 3.2 



We divide the proof into Part I and Part II. The statements in Part I are valid 



not only for C„ but for any § 



Apart from the separation into two solitons 

in which the 



in the final state, this already proves Theorem 3.2 for those g 
decomposition of B t 



Bk is multiplicity-free. Among the list of g n in question, 



Since the C„ case has the multiplicity, we 



such cases are A 2 n_i, B$ and Dn\ 
need Part II, which relies on the explicit result on the combinatorial R matrices 
for U' q (cP) in iHKOTfl . In the rest of Section | we shall write l^(a„)-crystals 
as B\ and C/g(g n _i)-crystals as B\ for distinction. 

Part I. For sufficiently large L, let ul — (L, 0, . . . ,0) S Bl be the highest 
weight element. Let 1Z : Bi <E> Bk — B k ® Bi be the combinatorial R matrix of 
C/g(§„_i)-crystals. The main claim in Part I is 

Proposition 3.8. For b\ € Bi,b2 G B k (l > k) and L\,L<z G Z>i, suppose 
there exist to £ Z>i, C\ € Bi and c 2 £ Bk such that 

n«>i2 s r-r\® L 3 



«Ufc(c 2 )«>|lJ ®t,(ci)®[lj 



(3.2) 



/io/rfs /or any t > to for some L' l7 L' 2 , L3 and L' 3 under the isomorphism Bf* ® 
Bi (g> • • • ® Bi ~ Bi <g> ■ ■ ■ <g> Bi <g) Bf*. De/me S : Bi ® B k B k ® Bi through 
this relation by S(b\ ® 62) = c 2 ® C\. Then lZ{b\ ® 62) = S(bi ® 62) implies 
Tt{fi{bi®b 2 )) = S{fi{bi®b 2 )) for each i = 1, . . . , n— 1. Similarly, €i(bi®b 2 ) = 
implies eiS(b x ® 62) = /or eac/i i = 1, . . . , n — 1. 



Lemma 3.9. For eac/i z = 1, . . 



, n — 1, we Ziaue a commutative diagram: 



e i+ i 



(Bi ® S fc ) U {0} -^!l> U {0}. 

XTie same relation holds also between fi and fi+i- 
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This is a simple corollary of Lemma |3.7| . Similarly we have 

Lemma 3.10. Let a%, ... , ar, G B\, i\, . . . ,i m be the subsequence of 1, ... ,L 
satisfying Pj = |T~| G Bi «/ j {ii,...,i m } (m < L), and pk = a lk (k 
1 m). TTien /or any t, t' G Z>o and £ G Z>i, 

1. /»+i(pi ® • • • ® p m ) = Pi ® ■ ■ ■ ® p' m 

=> /»+i(wf* ® ai <8 • • • <S> ai ® uf *') = uf * ® a' x <g> • • • <g a' L ® uf *', 

where a'={ P X * 3 = 

I 111 otherwise, 

2. fi+iiuf 1 ®a 1 ®---®a L ® uf *') = uf * ® a£ ® • • • <g> a' L <g> uf *' 
=>• fi+l(pi ® ■ • -®Pm) = Pi 8> • ■ • <S>p' m , where p' k = a' lk , 

3. f l+ i{uf ® ai ®---®a L ®uf') = O /i+i(pi®---®Pm) =0, 
/or eac/i i = 1, . . . ,n — 1. The same is true also for §j+i. 



Proof of Proposition jff.ffi Suppose TZ(b\ ® 6 2 ) = <S(bi ® 6 2 ) = c 2 ® ci. Thus 
(c 2 , Cx) G -Bfc x B; is related to (b\, 6 2 ) via (|3.2|). Set &i®&2 = /i(&i®&2) for some 
i = 1, . . . ,n— 1 and put c^tgic^ = 5(&i<g>& 2 ). We are to show 72.(6^ ®o 2 ) = c^<8)c^. 
From Lemma 3.9 one has 

(u(6i) ® »*(6a)) - tj(&i) ® ^( & 2 )- (3-3) 

00^1 I 1®^2 I 1®^3 

8) Zf (bi ) ®| 1 1 ®Zfc(o 2 )®[lJ ' and 

Pi ® • • ■ ®p m = k(bi) ® «fe(6 2 ) in Lemma |3.10| , one has 

/ i+ i (wf®[Tf Ll <g>zi(oi) ®|T| 8L2 ®? fe (o 2 )®[7f L3 ) 

= uf ® ® !,(&£) ® {Tf L2 ® z fc (6 2 ) ® 0^ 3 , 

where ( |3.3| ) is used. By sending ml to the right as in (|]^), this is equivalent to 

fi+i ® z fc (c 2 ) ®0^ 2 ®»i(ci) ®0® L3 ® uf 



= [Tf il ® « fe (4) ®[7f ia ®n(c[ ) ®[7f L3 ® u 



in terms of the (c 2 , c^) specified above. With the help of Lemma |3.1C| and |3.9| we 
can now go backwards to see that this implies /i+i(»fe(c 2 )(g)i/(ci)) = «fc(c 2 )®i;(cjj 
hence /i(c 2 ® ci) = c' 2 ® c^. Thus we have 1Z{b\ ® b' 2 ) = TZ(fi(b\ ® 6 2 )) = 
/i7?.(oi®6 2 ) = /i(c 2 (8>ci) = c^®^. It is very similar to verify that ei(6i®o 2 ) = 
implies e-iS(bi ® o 2 ) = for any i = 1, . . . , n — 1. □ 

The above proof is also valid in case [HHIKTT]. Viewed as U q (g n -x)- 
crystals with Kashiwara operators fi, e, (1 < i < n — 1), both crystals £?/ ® -Bfc 
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and Bk ® Bi decompose into connected components. Note that ( p.2[ ) obviously 
tells that [/q(g„_i)-weights of 61(8)62 and 5(61062) = c%®c\ are equal. Therefore 
apart from the separation into two solitons in the final state, Proposition |3 . 8| 
reduces the proof of Theorem [T^ to showing TZ{b\ ® 6 2 ) = S{b% ® 62) only 
for the ?7 g (0 n _i)-highest weight elements 61 ® 6 2 . In particular, if the tensor 
product decomposition of Bi ® Bk is multiplicity-free, it only remains to check 
the separation. 

Part II. Here we concentrate on the U' q (Cn^) automaton and prove the 
separation and 1Z(b\ ® 62) = <S(6i ® 62) directly for [/g(C n _i)-highest weight 
elements 61 ® 62 (n > 3). In the notation (x\, . . . , x n _x, S n _i, . . . ,5ix) € B). or 
Bi, they have the form 

{/,0 0)8(d,c I 1 ...,0 1 6)6B|8B tl / > 6 + c. (3.4) 

In what follows we always assume ^ > k and use the non-negative integers e and 
a defined by I — f + 2e and k = 2a + b + c + d. 



Proposition 3.11 ([HKOT|). Under the isomorphism TZ : Bi®B^ — B^®Bi 
of C^(C„_i) -crystals, the image of ( 3.4) is given by 

(k - 2e, 0, . . . , 0) <g> (d + I - k - y, c, 0, . . . , 0, 6 - y), y = mm(l -k,(b- d) + ), 

(3.5) 

ifa>e. In case a < e, it is given by 

{k — 2a, 0, . . . , 0) (g> (d + I — k — e + a — z, c, 0, . . . , 0, 6 + e — a — z), 
z = min(6 — d + e — a, I — k — e + a), 

if I — k>e~a>d— b, and 

(k - 2e + 2w, 0, . . . , 0) ® (d + I - k - w, c, 0, . . . , 0, 6 + w), 
w = mm(l — k, (2e — 2a — d + 6)+), 

otherwise. 



See Section 2.1 for the definition of the symbol (x)+. Below we employ the 
notation 

[s, t, u] — (x\ = L — s — t — 2u, X2 = s, x\ = t, other Xi, x~i = 0) £ -Bl, 

and always assume L 3> s,t,u. To save the space, the tensor product of 
U'(Cn )-crystals such as 

[s, t, uf j ® ^ P ® ® [a", t", u"] ® [if" ® [*', t', 
will simply be denoted by 

[s,t,u] j l p 2 q [s",t",u"] 3 r [s',t',u'] j '. 



From the results in [HKOT] we further derive two lemmas given below. 
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Lemma 3.12. If L p, g, s, t, u, we have 
[0,0,0] P2«~P+« [q,0,p], 

'l[s,t+l,u- 1] 
![«,*- 1,0] 
2 [s- 1,0,0] 
.1 [0,0,0] 



[s,t,u] 1 



u > 0, 

u = 0,t > 0, 
u = i = 0, s > 0, 
u = t = s = 0. 



Lemma 3.13. If L s, i, u, a, 6, c, d, we /iaue 
[s,t,u] l a 2 fc 3 C 2 d l 5 

^ J a c^s+t — a -^u—t— s-\-2a-\-b-\-c-\-d Ja+u 2^ — a +k ^ c 2^ [0 0] 

^ J a — a +2b+c -j^— t+2a— fe+d! Ja+u 2^~ o+fe ^s — lb—c-^rd |q q qi 



l a 2 



^a+ii ^s-\-t-\-u- 



-b—c—d c\2b-\-c-\-d—s nc 



¥2 d [0,0,0] 



J a 2^+ c +^ j^a+it 2^ 3 C 2 S_ c+d rQ q qi 

2 s ^u— s— t+2a+6+c+d ja+ti 2^ ^ c 2^ TQ 0] 



> 2 d 3 c 2S +t-a-2b-c+d j 0)0)() j 



where the conditions (I) -(VIII) are given by 



if (I), 
if (II), 

cm), 
*/(iv), 
*/(v), 

(VI), 



tf(vni), 



(I) 


t > 


a, 


s<2b + 


c, s + t<a + bH 


-c-f 


d, 




(II) 


t > 


a, 


s>2b + 


c. t < a — b + d, 








(III) 


t > 


a, 


s <2b + 


c, s + t > a + b-\ 


-c + 


■d, 




(IV) 


t > 


a, 


s>2b + 


c, t > a — b + d, 








(V) 


t < 


a. 


s + t < min(a + 2b + c, a - 




- c - 


\-d), 


(VI) 


t < 


a, 


b<d, s 


+ t > a + 26 + c, 








(VII) 


t < 


a, 


b > d, a 


+b+c+d< s+ 


< < 


a -+ 


-2b + c 


(VIII) 


t < 


a, 


b > d, s 


+ t > a + 2b + c. 









On the boundaries, the formulae in the adjacent domains are equal. The domain 
(II) is absent if b > d. S = a + b + c + d + t + 2u for (I, III), a-b + d+s + t + 2u 
for (II, IV, VI, VIII) , 2a + b + c + d + 2u for ( V, VII) . 



Proof of Theorem ff.i . Let b\ ® 62 be the highest weight el emen t ( |3.4[ ) 
and C2 ® c\ = 1Z(b\ ® 62) be its image specified in Proposition |3.1l[ Two 
solitons labeled by b\ and 62 emerge in our automaton as the patterns ii(b±) = 
l e 2* l e and i k (b 2 ) = l a 2 b 3 c 2 d l a , respectively. Thus S(h ® 6 2 ) can be found by 
computing the image of 



[0, 0, 0} j l e 2 f l e+m l a 2 b 3 C 2 d 11 • • • 



(3.6) 
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under the isomorphism B® 3 (8 {B± (g) _Bi 



• ) ~ (Bi <g) fl x <8) 



B 



for an interval m. This calculation can be done by using Lemma 3.12 and 
3.13 only. It branches into numerous sectors depending on the seven param- 
eters a, 6, c, d, f, e and m. (We use the dependent variables I — / + 2e and 
k = 2a + b + c + d simultaneously.) We have checked case by case that the 
collision always 'completes' within two time steps (j = 2) ending with the result 
l mi »fe(c2)l m2 «z(ci)ll ■ ■ ■ [0, 0, Op for some mi and m 2 . This establishes K = S 
hence Theorem [O. Below we illustrate a typical such calculation in the sector 



a > 2e, b > d, 



e<m</ + e — a — 2b — c. 



(3.7) 



as 



In this case c% ® C\ = 72.(t>i ® 62) is determined by Proposition 3.11 

ci = (I- k-b+2d,c,Q,... ,Q,d) e Bj, c 2 = (fc - 2e, 0, . . . , 0) e B fc . 
Thus we seek the patterns 



i k {c 2 ) = l e 2 fc 



k-2e 



»j(ci) = l c 



1° 



(3.8) 



to show up after a collision and separate from each other. Starting from ( p.6| 
with j = 3, one can rewrite it by means of Lemma 3.12 and (3.7) as 

[0, 0, 0] 2 l 4l l e 2 m - e [si, 0, 0]l a 2 fc 3 c 2 d ll 



with s\ = f + e — m and 81 = f + 2e. From ( p.7| ) one can apply Lemma 3.13 
(VIII) to transform this into 

[0, 0, Off 51 l e 2^ l a l° 2 2 d 3 c 2 d2 ll ... [0, 0, 0], 

where s' 2 = m — e + a + b + c + d,a 2 = a + b — d and d 2 = s\—a — 2b — c + d. This 
completes the calculation of the one time step in the automaton. The second 
time step can be worked out similarly. Lemma |3. 12 with (3.7) transforms the 
above into 



[0, 0, 0] l 5a l e 2 a - 2e [s 2 , 0, 0] l a2 2 d 3 c 2 d2 11 ... [0, 0, 0] 



where s % = s ' 2 — a + 2e = m + b + c + d + e and 8 2 = Si 
Lemma 3.13| (VI) is applied, leading to 



2e + So- This time 



[0 0] l^ 2 l e 2^ _ ^ e l^ 3 ^o,-\-b—d2^d^C2^l — k — b-\-2d 



ll...[0,0,0] 2 



(3.9) 



where 83 — a 2 — d + d 2 = / + e — m — b — c — d. This already contains the 
patterns ( |3.8D . Let us calculate one more time step to confirm that they are 
separating hereafter. Since 63 > 2e, (|3.9|) is isomorphic to 



l s *l e [k - 2e, 0, o]i^-^i a+b - d 2 d 3 c 2 l - k - b+2d 11 . . . [0, 0, 0] 2 , (3.10) 
where #4 = 82 + k. If k < 83 this becomes 

-|^4 ^e^k — 2e ^£3+/— k Ja+6— d^d^c^l — k — b+2d ^ |Q Q Qj 3 ^3 
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Compared with (|3.9|), the distance of the two patterns here has increased by 



I — k in accordance with the velocity of solitons in Theorem 3.1. In the other 
case k > S 3 one rewrites ( |3.10 ) as 



for which Lemma 3.13| (V) can be applied. The result coincides with (3.11). 



□ 



4 Discussion 

In this paper we have proposed a class of cellular automata associated with 
U' q (g n ). They are essentially solvable vertex models at q = in the vicinity of 
the ferromagnetic vacuum. They exhibit soliton behavior stated in Theorem 



ll and Conjecture 3.3 



Behavior of tensor products of crystals in the vicinity of ferromagnetic vac- 
uum has not been investigated in detail so far. To understand it better will be a 



key to prove Conjecture 3.3 and to clarify the relevance to the subalgebra g n -i, 
which has also been observed in the RSOS models in the ferromagnetic-like 
'regime IF pR|,|i^|. 



The automata considered in this paper are associated with an / — > 00 limit 
R' : Bu (g> Bi — ► Bx (g> B^ of the ordinary combinatorial R matrices R : Bi® Bi ~ 
Bx® B[. To vary Bi from site to site and to replace R' with the corresponding 



combinatorial R matrices is a natural generalization. See [ HHIKTT ] . 

We close by raising a few questions. What kind of soliton equations can pos- 
sibly be related to the U' q (g n ) automata in general? Is it possible to derive them 
conceptually from the associated vertex models, bypassing the route; vertex 

, , q— >0 tropical variable change . . , , . n T . , rn • 

models ► automata > soliton equations ! Is it efficient to 

seek a piecewise analytic form of the combinatorial R matrices Bi®Bk ~ B)~<2)Bi 
by ultra-discretizing the relevant solitons solutions? Is it possible to extract in- 
formations on the associated energy function from the automata? 

A The map B! : B^® B 1 -> B 1 (g) B^ 

Let gi : B{ — > B^ be the embedding of Bi as a set as mentioned in Section ^. By 
the definition, for any u £ B^ there exists the unique element ui £ Bi such that 
giiui) = u for each I which is sufficiently large. Consider the composition: 

Bl, x B x 91 1 ® ld : Bi ® Bx {Bi ® Bx) U {0} x Bx) U {0} 

(u,b) 1 — > ui®b 1 — > u'^b 1 1 — > (jgi{u'i),V), 
where of course either u[ — ui or b' — b, and we interpret (<?; x id)(0) = 0. For 
all the Bi considered in this paper it is easy to sec that this is independent of I 
if I is sufficiently large. We let B^ ® Bx denote the set B^ x Bx equipped with 
the actions e-, /• determined as above with sufficiently large I. We shall simply 
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write it as e-, f[ : B^ (g> B\ — > (B^ <g> Bi) U {0}. Similarly we define Bi B^ and 
ej, # : fli ® flu-* (Bi ® B h ) U {0}. 

Given u G B^ and G Bi, let c = c(l,u,b) G Bi and u = v(l,u,b) G B; be 
the elements determined by 

under the combinatorial B matrix B : B; <g> Bi ~ Bi ® B;. Then we have 

Conjecture A.l. For any fixed u G B^ andb G Bi, </ie elements c{l,u,b) G Bi 
and gi{y{l,u,b)^ G B^ are independent of I fori sufficiently large. 

Assuming the conjecture we define the map R! : B[, <g> B\ — > Bi ® B[, to be 
the composition: 

R : B^® B 1 9 ' 1(S ' d - Bi ® B x ~ Bx <g> B; Bx ® B^ 

m <g> 6 i — > ui®b ~ c®u i — ► c®gi(v), 
with Z sufficiently large. Obviously this is invertible. Moreover it commutes 
with and although this fact is not used in the main text. For this can 
be seen from the commutative diagram (f- case is completely parallel.) 

B ^ Bl B®Bx — «-> Bi ® B; B®B t 



<1 "i 



B^0Bi g ' B;®B! — £-> Bi ig) B; _==S2»* B (g) Bt,, 

where the left and the right squares are the definitions of e£ (for elements not 
annihilated by e-). The map B' indeed has the property (I) in Section [|. The 
property (II) is also valid for and Cn ■ We conjecture it for all the cases 
considered in this paper. 



Note: While writing the paper the authors learned from FOY that the energy 



of combinatorial R matrices for U' q (An ) is encoded in the phase shift of soliton 
scattering in the automaton. We thank Yasuhiko Yamada for communicating 
their result. 



References 

[B] R.J. Baxter, Exactly solved models in statistical mechanics, Academic Press, 
London (1982). 

[BR] V. V. Bazhanov and N. Yu. Reshctikhin, Restricted solid- on- solid models 
connected with simply laced algebras and conformal field theory, J. Phys. A. 
Math. Gen. 23 (1990) 1477-1492. 

[DJKMO] E. Date, M. Jimbo, A. Kuniba, T. Miwa and M. Okado, One- 
dimensional configuration sums in vertex models and affvne Lie algebra char- 
acters, Lett. Math. Phys. 17 (1989) 69-77. 



27 



[DJM] E. Date, M. Jimbo and T. Miwa, Representations of U q (Ql(n, C)) at 
q = and the Robins on- Schensted correspondence, "Memorial Volume for 
Vadim Kniznik, Physics and Mathematics of Strings" , eds. L. Brink, D. 
Friedan and A. M. Polyakov, World Scientific, Singapore (1990) 185-211. 

[FPS] A. S. Fokas, E. P. Papadopoulou and Y. G. Saridakis, Soliton cellular 
automata, Physica D41 (1990) 287-321. 

[FOY] Fukuda, M. Okado, Y. Yamada, private communication. 

[HHIKTT] G. Hatayama, K. Hikami, R. Inoue, A. Kuniba, T. Takagi and 
T. Tokihiro, The Automata related to crystals of symmetric tensors, 
preprint, |math.QA/9912209| . 

[HKKOT] G. Hatayama, Y. Koga, A. Kuniba, M. Okado and T. Takagi, Finite 



crystals and paths, preprint, math.QA/9901082 



[HKOT] G. Hatayama, A. Kuniba, M. Okado and T. Takagi, Combinato- 
rial R matrices for a family of crystals: Cn^ and cases, preprint, 



math.QA/9909068 



[HI] K. Hikami and R. Inoue, private communication. 

[HIK] K. Hikami, R. Inoue and Y. Komori, Crystallization of the Bogoyavlensky 
lattice, J. Phys. Soc. Jpn. 68 (1999) 2234-2240. 

[Kas] M. Kashiwara, Crystalizing the q-analogue of universal enveloping alge- 
bras, Commun. Math. Phys. 133 (1990) 249-260. 

[KKM] S-J. Kang, M. Kashiwara and K. C. Misra, Crystal bases of Verma 
modules for quantum affine Lie algebras, Compositio Math. 92 (1994) 299- 
325. 

[KMN1] S-J. Kang, M. Kashiwara, K. C. Misra, T. Miwa, T. Nakashima and 
A. Nakayashiki, Affine crystals and vertex models, Int. J. Mod. Phys. A 7 
(suppl. 1A), (1992) 449-484. 

[KMN2] S-J. Kang, M. Kashiwara, K. C. Misra, T. Miwa, T. Nakashima and 
A. Nakayashiki, Perfect crystals of quantum affine Lie algebras, Duke Math. 
J. 68 (1992) 499-607. 

[KN] M. Kashiwara and T. Nakashima, Crystal graphs for representations of 
the q-analogue of classical Lie algebras, RIMS preprint (1991) to appear in 
J. Alg. 

[Ki] A. N. Kirillov, Discrete Hirota equations, Schutzenberger involution and 
rational representations of extended affine symmetric group, talk at "Devel- 
opments in Infinite Analysis II", Okinawa, December 1998. 

[KS] P. P. Kulish and E. K. Sklyanin, Quantum spectral transform method. 
Recent developments, Lect. Note. Phys. 151 Springer (1982) 61-119. 



28 



[K] A. Kuniba, Thermodynamics of U q (X^) Bethe ansatz system with q a root 
of unity, Nucl. Phys. B389 (1993) 209-244. 

[MSTTT] J. Matsukidaira, J. Satsuma, D. Takahashi, T. Tokihiro and M. Torii, 
Toda-type cellular automaton and its N-soliton solution, Phys. Lett. A 225 
(1997) 287-295. 

[N] T. Nakashima, Crystal base and a generalization of the Littlewood- Richard- 
son rule for the classical Lie algebras, Commun. Math. Phys. 154 (1993) 
215-243. 

[NY] A. Nakayashiki and Y. Yamada, Kostka polynomials and energy functions 
in solvable lattice models, Selecta Mathematica, New Ser. 3 (1997) 547-599. 

[PF] T.S. Papathcodorou and A. S. Fokas, Evolution theory, periodic particles, 
and solitons in cellular automata, Stud. Appl. Math. 80 (1989) 165-182. 

[PAS] T.S. Papatheodorou, M. J. Ablowitz and Y. G. Saridakis, A rule for fast 
computation and analysis of soliton automata, Stud. Appl. Math. 79 (1988) 
173-184. 

[PST] J. K. Park, K. Steiglitz and W. P. Thurston, Soliton-like behavior in 
automata, Physica D19 (1986) 423-432. 

[STF] E. K. Sklyanin, L. A. Takhatajan, L. D. Faddecv, Quantum inverse prob- 
lem method L, Theor. Math. Phys. 40 (1980) 688-706. 

[TS] D. Takahashi and J. Satsuma, A soliton cellular automaton, J. Phys. Soc. 
Jpn. 59 (1990) 3514-3519. 

[T] D. Takahashi, One some soliton systems defined by using boxes and balls, 
Proceedings of the International Symposium on Nonlinear Theory and Its 
Applications (NOLTA '93), (1993) 555-558. 

[TTMS] T. Tokihiro, D. Takahashi, J. Matsukidaira and J. Satsuma, From soli- 
ton equations to integrable cellular automata through a limitting procedure, 
Phys. Rev. Lett. 76 (1996) 3247-3250. 

[TNS] T. Tokihiro, A. Nagai and J. Satsuma, Proof of solitonical nature of box 
and ball systems by means of inverse ultra-discretization, preprint. 

[W] S. Wolfram, ed. Theory and applications of cellular automata, World Sci- 
entific, Singapore (1986). 



29 



